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The NASA Double Asteroid Redirection Test (DART) mission performed a kinetic impact 
on asteroid Dimorphos, the satellite of the binary asteroid (65803) Didymos, at 23:14 UTC 
on September 26, 2022 as a planetary defense test’. DART was the first hypervelocity 
impact experiment on an asteroid at size and velocity scales relevant to planetary 
defense, intended to validate kinetic impact as a means of asteroid deflection. Here we 
report the first determination of the momentum transferred to an asteroid by kinetic 
impact. Based on the change in the binary orbit period’, we find an instantaneous 
reduction in Dimorphos's along-track orbital velocity component of 2.70 + 0.10 mm s", 
indicating enhanced momentum transfer due to recoil from ejecta streams produced by 
the impact?^^, For a Dimorphos bulk density range of 1,500 to 3,300 kg m^, we find that 
the expected value of the momentum enhancement factor, 5, ranges between 2.2 and 4.9, 
depending on the mass of Dimorphos. If Dimorphos and Didymos are assumed to have 
equal densities of 2,400 kg m^,  — 3.61*0123 (10). These f values indicate that 
significantly more momentum was transferred to Dimorphos from the escaping impact 
ejecta than was incident with DART. Therefore, the DART kinetic impact was highly 


effective in deflecting the asteroid Dimorphos. 


Observations from the DART spacecraft on approach found Dimorphos to be an oblate spheroid 
with a boulder-strewn surface, and the spacecraft impacted within 25 m of the center-of-figure'. 
Ejecta from the DART impact were observed in situ by the Italian Space Agency's Light Italian 
Cubesat for Imaging of Asteroids (LICIACube) spacecraft, which performed a flyby of 
Dimorphos with closest approach about 168 s after the DART impact?. The impact ejecta were 
further observed by Earth- and space-based telescopes, revealing ejecta streams and dust tails 
similar to those seen in active asteroids thought to be triggered by natural impacts?9. Ground- 
based telescopes and radar determined that the DART impact reduced the binary orbit period 
by 33.0 + 1.0 (3e) minutes?. 


As a planetary defense test mission, a key objective of DART is to determine the amount of 


momentum transferred to the target body relative to the incident momentum of the spacecraft, 


quantified by the momentum enhancement factor £ (e.g., refs. 478), which is defined by the 


momentum balance of the kinetic impact, 
MA = mU + m(B-1)(E-U)E. (1) 


Here, M is the mass of Dimorphos, A7 is the impact-induced change in Dimorphos's orbital 
velocity, m is DART's mass at impact, U is DART's velocity relative to Dimorphos at impact, and 
E is the net ejecta momentum direction. MAÐ is the momentum transferred to Dimorphos, mU is 
DART's incident momentum, and the final term in the equation is the ejecta's net momentum 
written in terms of the spacecraft incident momentum. In this formulation, f is the ratio of actual 
imparted momentum to the impactor's momentum in the direction of the net ejecta momentum. 
Although previous works have defined f using the impactor's momentum in the surface normal 
direction??, we elect to instead use the ejecta direction as our reference in order for the result to 
be independent of the surface topography. These definitions are equivalent in the case where 
the ejecta direction is in the surface normal direction. A B near 1 would suggest that ejecta recoil 
had made only a negligible contribution to the momentum transfer. A B > 2 would mean that the 


ejecta momentum contribution exceeded the incident momentum from DART. 


The full Av cannot be determined with the available information'?, but its component along 
Dimorphos's orbital velocity direction, referred to as the along-track direction, can be estimated 
from available data including Dimorphos's orbit period change. To express f in terms of the 
along-track component of AV, we take the scalar product of (1) with the unit vector êr in the 
along-track direction. Solving for p} yields, 

M — 

z (Avr) - (U-ér) 

B m (£-U)(E-er) (2) 

For the remainder of this work, we refer to the along-track component of Dimorphos's velocity 
change, Av - êr, as Avy. Fig. 1 shows the geometry of the DART impact, including the nominal 
ellipsoidal shapes used for Didymos and Dimorphos in our analysis, and the nominal 


orientations of U, êr, and E at the time of impact. 


The major unknowns in calculating £ are Avy, M, and E. We first employ a Monte Carlo 
approach to produce a distribution for Avy consistent with the measured period change that 
incorporates the various uncertainties involved. We sample many possible combinations of 


Didymos system parameters, including the ellipsoid shape extents of the asteroids, pre-impact 


orbit separation distance between the two asteroids’ centers of mass (i.e., Dimorphos’s pre- 
impact orbit radius), pre- and post-impact orbit periods, and net ejecta momentum direction ÊF. 
We use the full two-body problem code (GUBAS"', see Methods) which implements coupled 
rotational and orbital dynamics to numerically determine Av; for each sampled combination of 
input parameters. Coupled dynamics are necessary because of the non-spherical shapes of 
Didymos and Dimorphos and their close proximity relative to their sizes. A range of values for M 
is generated by combining the volumes of the sampled ellipsoid shape parameters with values 
for Dimorphos’s density. Since Dimorphos’s density has not been directly measured and has a 
large uncertainty, we treat it as an independent variable and uniformly sample a wide range of 
possible values between 1,500 and 3,300 kg m^, a range that encompasses the 30 
uncertainty’. Using a technique modified from that of ref. ^ (see Methods), we apply 
observations of the ejecta via Hubble and LICIACube data to obtain a preliminary measurement 
of the axis of the ejecta cone geometry. The cone axis direction is identical to E assuming the 
ejecta plume holds the momentum uniformly, and we find E points toward a right ascension 
(RA) and declination (Dec) of 138° and +13°, respectively (see Extended Data Fig. 1). We 
assign a conservative uncertainty of 15° around this direction. Finally, £ also depends on 
DART's mass and impact velocity, as well as Didymos's pole orientation’. Those quantities have 
negligibly small uncertainties relative to those of the other parameters discussed previously and 
are therefore treated as fixed values (not sampled). See Methods for additional details on the 
Monte Carlo analysis, and Extended Data Table 1 for a list of parameters and uncertainties, and 


Extended Data Table 2 for the covariances that were utilized. 


We find that Av; = —2.70 + 0.10 (1c) mm s^, based on the observed impact-induced period 
change of —33.0 + 1.0 (30) minutes and the shapes and separation of Didymos and 
Dimorphos'?. Fig. 2 shows the distribution of Av; values from the Monte Carlo analysis, along 
with the fitted mean and standard deviation. The resulting spread of 6 values as a function of 
Dimorphos's density, calculated via Eq. (2), is presented in Fig. 3, along with linear fits for the 
mean f versus density trend and its 10 confidence intervals. The linear-fit slope is expressed as 
a scale factor on the ratio of density to the nominal value of 2,400 kg m? (ref. 1). For that 
nominal Dimorphos density, at which Dimorphos and Didymos would have approximately equal 
densities, £ = 3.61*022 with 1e confidence. The mean £ ranges between 2.2-4.9 as a function 
of density across the range of 1,500 to 3,300 kg m^, and overall, 6 ranges between 1.9—5.5 


with 3e confidence. 


Our result for £ is consistent with numerical simulations'?-? and laboratory experiments??? of 
kinetic impacts, which have consistently indicated that 6 is expected to fall between about 1 and 
6. However, non-unique combinations of asteroid mechanical properties (e.g., cohesive 
strength, porosity, and friction angle) can produce similar values of f in impact simulations '?. 
Future studies that combine estimates of 6 with additional constraints from the DART impact 
site geology! and ejecta observations?? will provide greater insight into Dimorphos's material 
properties. In addition, ESA's Hera mission? is planned to arrive at the Didymos system in late 
2026. By measuring Dimorphos's mass and other orbital properties, Hera will allow us to 


significantly improve the accuracy and precision of the p determination. 


DART's impact demonstrates that the momentum transfer to a target asteroid can significantly 
exceed the incident momentum of the kinetic impactor, validating the effectiveness of kinetic 
impact for preventing future asteroid strikes on the Earth. The value of $} from a kinetic impact is 
key to informing the strategy of a kinetic impactor mission (or missions) to mitigate a future 
asteroid impact threat to Earth??. Should £ turn out to be greater than 2 across a wide range of 
asteroid types, it would mean significant performance improvements for kinetic impactor 
asteroid deflection missions. If f > 2, as opposed to } ~ 1, then the same sized kinetic impactor 
could deflect a given asteroid with less warning time, or deflect a larger asteroid with a given 


warning time than it could otherwise. 
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Figure 1. Schematic of the DART impact geometry on Dimorphos. The pre-impact orbit is 
shown with a solid line around Didymos. The dashed line sketches the orbit change due to the 


impact. Orbits are drawn approximately to scale. The positive pole direction of Didymos is h 


(pointing down in the bottom panel). DART's incident direction is U, the net ejecta momentum 
direction is E (which points to a right ascension (RA) and declination (Dec) of 138° and +13°, 
respectively), and the direction of Dimorphos's orbital motion, referred to as the along-track 
direction, is êr. The relative positions of the Sun and the Earth are also indicated. The upper 
panel shows the view from Didymos's negative pole direction, while the lower panel provides a 


perspective view. 
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Avr = —2.70 £0.10(1c) mm s^! 
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Figure 2. Probability distribution of Avy, the along-track component of the change in 
Dimorphos's velocity induced by DART's impact, generated by our Monte Carlo analysis that 
samples over input parameter uncertainties. The histogram consists of 100,000 Monte Carlo 


samples and is normalized to an area of unity. A Gaussian fit to the distribution indicates a 


mean Av, of —2.70 mm s^! with a standard deviation of 0.10 mm s^. 
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Figure 3. £ as a function of Dimorphos's bulk density pg, from the dynamical Monte Carlo 
analysis. Individual samples are plotted as points, while the linear fit for the mean f is plotted as 
the solid line and the dotted lines show the 10 confidence interval. The color bar indicates the 
mass of Dimorphos corresponding to each Monte Carlo sample, which is determined by bulk 
density and the volume. The density range shown corresponds to the 3e range of the Didymos 
system density, while the shaded region highlights the 16 range’. If the density of Dimorphos 
were 2,400 kg m^, the densities of Didymos and Dimorphos would be the same as the system 
density, and f = 3.61*912 (10). For context, the densities of 3 other S-type near-Earth asteroids 
are in the range shown: 433 Eros?? at 2,670 + 30 kg m^5; 25143 Itokawa?? at 1,900 + 130 kg 
m~; 66391 Moshup?? at 1,970 + 240 kg m. 


12 


Methods 


Numerical Determination of 4v; and £ 

Several parameters affect the value of f as presented in Eq. (2): Avr, M, and E. The along- 
track velocity change, Avr depends on orbit period change, pre-impact semimajor axis, and the 
shapes of Didymos and Dimorphos, while M depends on Dimorphos's shape and bulk density 
(which was not measured). £ is the only parameter that is directly observed, but it still has 
considerable uncertainty. Thus, there are 12 total unknown input parameters: 3 axis lengths for 
Didymos's ellipsoidal shape (4x, Ay, Az), 3 axial lengths for Dimorphos's ellipsoidal shape 
(Bx, By, Bz), Dimorphos's bulk density pg, the pre-impact orbit semimajor axis apre, pre-impact 
orbit period P,e, and post-impact orbit period P,,,., and 2 angles to define the ejecta 
momentum direction vector (E). Extended Data Table 1 lists these input parameter values and 
their uncertainties, along with additional known quantities needed to calculate 6. To account for 
this large set of input uncertainties, we use a Monte Carlo approach in which 100,000 possible 
cases are generated by randomly sampling the input parameters within their uncertainties. We 
assume the DART spacecraft mass, DART impact velocity vector, and Dimorphos's orbital 
velocity direction (referred to as the along-track direction) are all known precisely because their 


uncertainties are negligibly small compared to the uncertainties of the other input parameters. 


The pre-impact orbit semimajor axis, pre-impact orbit period, and post-impact orbit period are 
sampled as a multivariate Gaussian distribution using the mean values and covariance matrix 
from the “N22+” solution (refs. 3132 of ref. 2; see Extended Data Tables 1 and 2). This accounts 
for the small correlations between those three parameters. The physical extents of Didymos and 
Dimorphos from ref. ' are sampled uniformly, as those uncertainties are not Gaussian (see 
Extended Data Table 1). £ depends strongly on Dimorphos's mass, but the mass is poorly 
constrained because Dimorphos's bulk density has not been directly measured'. Therefore, we 
treat density as the independent variable, sample it uniformly, and report £ as a function of 


Dimorphos's density. 


For each Monte Carlo sample for the Didymos system, a secant search algorithm (a finite- 
difference Newton's Method) described in ref. °° is first used to compute the density of Didymos 
required to reproduce the sampled pre-impact orbit period, given the sampled pre-impact orbit 
semimajor axis, body shapes, and Dimorphos's density. Then, a second secant search 


algorithm is used to determine the Av; required to achieve the sampled post-impact orbit 
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period. We match the pre- and post-impact orbit periods because these are directly measured 
by ground-based observations and are thus the best-constrained parameters of the system?. 
Given the non-Keplerian nature of the Didymos system, we use the General Use Binary 
Asteroid Simulator (GUBAS) to numerically propagate the binary asteroid dynamics. GUBAS is 
a well-tested full two-body problem (F2BP) code that can model the mutual gravitational 
interactions between two arbitrarily shaped rigid bodies with uniform mass distributions'!!?^, 
GUBAS has been benchmarked against other F2BP codes% and used extensively in previous 
dynamical studies of the Didymos system (e.g., '°°3:°°). Finally, the mass of Dimorphos is 
calculated from its ellipsoidal shape and Dimorphos's density. This mass, along with the 
computed Av; and sampled net ejecta momentum direction, are provided as inputs to Eq. (2) to 
calculate the value of corresponding to each of the 100,000 realizations of the system. For a 
discussion on estimating £, see the Ejecta Plume Direction section below. The process 


described herein is summarized graphically in Extended Data Fig. 1. 


The convergence criteria on both secant algorithms are set such that the simulated orbit period 
matches the desired orbit period to an accuracy 10 times better than the uncertainty on the 
measurements themselves. The numerical simulations measure the average orbit period of 
Dimorphos in an inertial frame over 30 days to account for small fluctuations in the mutual orbit 
period resulting from spin-orbit coupling?9. Our selection of 100,000 as the number of samples 
to use in the Monte Carlo analysis was informed by calculating an estimated minimum 
necessary number of samples from the Central Limit Theorem and then testing sample sizes 


near that estimated value. The £ estimate results are well converged with 100,000 samples. 


In the numerical simulations, both Didymos and Dimorphos are modeled as triaxial ellipsoids 
with physical extents from ref. 1. Images from DRACO and LICIACube showed that both 
Didymos and Dimorphos have an oblate spheroid shape’. There is no advantage to using more 
sophisticated shape models while the internal mass distributions of the bodies are unknown. 
Instead, the ellipsoidal approximation allows for easy sampling of a range of plausible moments 
of inertia as a proxy for different internal density distributions. For example, given the current 
uncertainties in Didymos's physical extents!, sampling over the given range of ellipsoidal 
shapes results in a range of plausible 2nd-order gravity terms (analogous to the spherical 
harmonic terms J2, Cz», etc.), which play an important role in the system's dynamics due to the 
tight separation of the binary components. Neglecting their shapes and assuming Keplerian 


dynamics results in Avr = —2.86 + 0.095(16), whereas GUBAS's 2nd-order gravity model finds 
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Avr = —2.70 + 0.10(10). Although 4th-order dynamics influence higher-order dynamical 
effects???9. we find that it comes with a significantly higher computational cost yet plays a 
negligible role in determining Avy. A smaller batch (due to increased computational cost) of 
~4,000 runs was conducted with 4th-order dynamics, which resulted in Avy = —2.68 + 0.10(10), 
suggesting the 2nd-order dynamics model is appropriate for determining Av; given the current 
uncertainties in the orbit solution and body shapes. This result was also independently verified 
using analytical models, accounting for Didymos's gravitational quadrupole, which agreed within 
a few percent of the 2nd-order numerical results, as expected given their dynamical 


approximations. 


We do not sample the rotation period of Dimorphos, as it is assumed to be equal to the pre- 
impact orbit period prior to the impact, with reasoning as follows. A measured inward orbit 
semimajor axis drift?" implies the system is evolving under the influence of the BYORP effect*®, 
which requires a secondary in near-synchronous rotation. Additionally, radar images constrain 
Dimorphos's spin period to be within 3 hours of the synchronous rate?!. Recent models for tidal 
dissipation in binary asteroids suggest that any free libration would dissipate on hundred-year 
timescales??, making any significant free libration unlikely given the timescales for excitation 
mechanisms such as close planetary encounters and natural impacts '?. Furthermore, 
Dimorphos's pre-impact eccentricity is constrained to be less than 0.033137, putting the 
maximum possible forced libration amplitude^? at ~0.5°. Although Dimorphos's rotation state is 
not precisely determined by DART, this body of evidence suggests that Dimorphos was likely in 


near-synchronous rotation and on a nearly circular orbit prior to the DART impact. 


Our model further assumes all momentum is transferred instantaneously, since earlier work 
showed the time duration of the momentum transfer has a negligible effect on the resulting 
dynamics’. The instantaneous torque on Dimorphos due to DART's slightly off-center impact is 
also neglected as the corresponding change in Dimorphos's rotation state is small compared to 
that arising from exciting Dimorphos's eccentricity and libration by the impact'??5^!, Finally, the 
effects of reshaping and mass loss due to cratering and ejecta are also neglected, as these 
effects are expected to be smaller in magnitude than the current ^1 min uncertainty on the post- 
impact orbit period^? and will remain poorly constrained until the Hera mission characterizes the 
Didymos system in 20277°. We leave these higher-order effects for future work once the post- 


impact orbit solution is refined further. 
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Ejecta plume direction 


We use observations of the ejecta plume to determine the ejecta momentum direction £. The 
conical ejecta plume was imaged by the LICIACube LUKE camera? and the Hubble Space 
Telescope (HST)*. We apply a technique used to derive cometary spin poles? to estimate the 
orientation of the ejecta cone axis. Although it is possible to have an asymmetric distribution of 
ejecta momentum (mass and velocity) within the cone, we assume the cone to be axially 
symmetric. The approach applies the ejecta cone's bright edges (if captured in an image) to 
compute the apparent direction of the cone axis projected onto the sky, which is assumed to be 


the middle of the edges. 


For a LICIACube observation, the projected cone axis defines the LICIACube-axis plane in 
inertial space that contains the line-of-sight and the projected axis. The cone axis can lie 
anywhere in this plane. The analogous plane HST-axis is defined from early HST images of the 
plume (those taken within 2 h after the impact) that show similar radial velocity as the ejecta in 
the LICIACube images, indicating it is likely the same ejecta material observed at a larger 
spatial scale. The intersection of these planes defines the cone axis orientation in three 
dimensions, but unfortunately the LICIACube-axis and HST-axis planes are nearly parallel. 
Thus, these observations do not provide a unique solution but they constrain the axis orientation 
to a narrow swath of the sky (see Extended Data Fig. 2). However, LICIACube LUKE images 
resolved the ejecta cone and the cone morphology over a large range of viewing angles during 
the flyby, further constraining the cone axis orientation?. During the approach to Dimorphos, the 
cone was pointed toward LICIACube, with the ejecta obscuring parts of Dimorphos. During 
recession from Dimorphos, the cone pointed away from Dimorphos and revealed it in silhouette 
(see Extended Data Fig. 3). The tightest constraint on cone orientation would have come from 
closest-approach images, with the cone axis in the plane-of-sky, as the cone transitioned from 
pointing toward the observer to pointing away. Unfortunately, Dimorphos and the ejecta cone 
were outside the LUKE field-of-view for 13 s around closest approach, and we lack images from 


the transition. 


The resolved LICIACube images are used to eliminate portions of the swath in Extended Data 
Fig. 2 where the observed cone morphology is inconsistent with those axis orientations. For 
example, half of the swath is rejected because the axis would be pointed in the opposite 
direction of what was observed. We also exclude orientations in which the cone would point too 


close to the line-of-sight during the approach or recession of LICIACube. We find the axis 
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orientation to be (RA, Dec) = (138°, +13°). We assign conservative uncertainties of ~15° in all 


directions based on the angular extents of region 5 in Extended Data Fig. 2. 
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Data availability 


The dynamical simulations were carried out using GUBAS, which is publicly available on Github 
(https://github.com/meyeralexj/gubas). The Dimorphos orbital velocity direction vector 
components presented in Extended Data Table 1 were computed using the 
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https://dart.jhuapl.edu/SPICE kernels/spk. The DART incident velocity vector components 
presented in Extended Data Table 1 were computed using those two files in combination with 
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availability at that URL is planned until summer 2023, after which those data may be found at 


https://naif.jpl.nasa.gov/naif/data.html. 
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Extended Data 


Extended Data Table 1: Values and uncertainties used for numerical simulations. All 
vectors are reported in the Earth Mean Equator J2000 (EME J2000) coordinate frame, at the 
impact time of September 26, 2022, 23:14:24.183 UTC". For Gaussian uncertainties we report 
the 10 uncertainties. For uniform uncertainties we report the median and the range of possible 


values. 


Quantity Value Uncertainty Assumed in 

Monte Carlo simulations 
DART mass atimpact,m | 579.4 Not considered Actual uncertainty! +0.7 
[kg] 
DART incident velocity 3.57399 Not considered See Data Availability for source 
vector x [km/s], Ux information. 
DART incident velocity -4.64106 | Not considered See Data Availability for source 
vector y [km/s], Uy information. 
DART incident velocity -1.85622 | Not considered See Data Availability for source 
vector z [km/s], U; information. 
Orbit period pre-impact, | 11.92148 | +0.000044 1o, Gaussian, reported by the “N22+” 
Fore [hr] solution in ref. 2 
Orbit period post-impact, | 11.372 +0.0055 1o, Gaussian, reported by the “N22+” 
Ppost [hr] solution in ref. 2 
Pre-impact orbit +0.035 1o, Gaussian, reported by the “N22+” 
separation (semi-major solution in ref. 2 
axis), Ayre [km] " : 


Didymos major axis, A, 


[m] 


Didymos intermediate 
axis, Ay [m] 


Didymos minor axis, A, 


[m] 


sampled uniformly 
sampled uniformly 


Dimorphos major axis, B, 


[m] 
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Dimorphos intermediate | 174 +4 sampled uniformly! 
axis, By [m] 

Dimorphos minor axis, B, | 116 +2 sampled uniformly! 
[m] 


Dimorphos density, pg 2400 +900 Sampled uniformly. Three times the 
-3 
[Kem] density uncertainty on the combined 
Didymos system from ref. ! 


Ejecta cone direction x, -0.72410 | Uncertainty is within a Sample uniformly throughout a cone 
15? cone around this of directions up to 15? away nominal 


Ejecta cone direction z, 0.22495 
E; 
Dimorphos orbital -0.689282 | Not considered See Data Availability for source 
velocity direction x, êr ; i 

X information. 
Dimorphos orbital 0.716645 | Not considered See Data Availability for source 
velocity direction y, êr J M Cum. 
Dimorphos orbital 0.106350 | Not considered See Data Availability for source 
velocity direction z, êr , i 

a information. 


0.65198 nominal vector RA,Dec = 138°,+13° 
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Extended Data Table 2: Covariance matrix. The covariances used to sample semimajor axis, 
apre, and pre- and post-impact orbit periods. The orbital solution from ref. ??? fits the pre-impact 
mean motion at the impact epoch, npre, and the change in mean motion due to the DART 
impact, An. Once these parameters are sampled, they are turned into a pre- and post-impact 
orbit period by the relation Pyre = 21: /n,,, and Ppost = 2m/(n,,* + An). The covariance matrix 


is constructed using units of km for apre and rad/s for npre and An. 


pow Lo Lo 


Secant Search Algorithm 


Iterate Avy to match Prost 
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Extended Data Figure 1. Outline of algorithm used to calculate P. The Monte Carlo variables 
are outlined by dashed lines and are defined as follows: pre-impact orbit period P,,.; pre-impact 
orbit semimajor axis apre; Didymos ellipsoid extents Ay, 4), Az; Dimorphos ellipsoid extents 

B, By, B;; Dimorphos density pg; post-impact orbit period Ppost; and net ejecta momentum 
direction £. First, a secant algorithm iterates Didymos density p4 to match Pyre- Next, another 
secant algorithm iterates the along-track change in Dimorphos's velocity Av; to match Ppost- 
Finally, M is calculated using the ellipsoid extents and density of Dimorphos, and then combined 


with Av; and Ê to calculate £. 


LICIACube 


Region @ 


Declination (deg) 


300 270 240 210 180 150 120 90 60 30 
Right Ascension (deg) 


Extended Data Figure 2. Ejecta cone orientation lies in the swaths of sky (black lines) defined 
by HST and LICIACube observations. The light-blue envelope outlines the axis position 
uncertainty in the direction measured in the sky plane. Red lines divide the along-plane swaths 
into regions that are excluded based on cone morphology in LICIACube images: 1) and 2) are 
excluded because the ejecta cone would point in the opposite direction from what is observed; 
3) is excluded because the axis would lie too close to the sky plane; 4) is excluded because the 
axis would lie too close to the line-of-sight; and 5) is the expected region for the axis orientation. 
The yellow dot denotes the best solution (RA,Dec) = [138° ,+13°] with the dark-blue envelope 


showing the extent of possible solutions. The red square is the direction of the incoming DART 
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trajectory [128°,+18°] and the green triangle shows Dimorphos's velocity vector [134°,+5°]. The 
LICIACube swath is defined for the +175 s image shown in Extended Data Fig. 3. 


Extended Data Figure 3. Two LICIACube LUKE images showing the ejecta morphology that 
were used to reduce the possible axis orientation solutions. The left panel shows an approach 
observation, 156 seconds after impact, with the ejecta in front of and partially obscuring 
Dimorphos. The right panel shows the ejecta morphology after close approach, 175 seconds 
after impact, with Dimorphos silhouetted against the ejecta cone. The images show the red 
channel from frames LICIACUBE LUKE L2 1664234220 00005 01 and 
LICIACUBE LUKE L2 1664234239 01003 01. The bright object in the upper corner of each 


image is Didymos. 
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